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Abstract 

A dynamical coupled channel approach is used to study ir and o;-meson production induced 
by pions and photons scattering from the proton. Six intermediate channels including vrA^, nN, 
7rA, <jN, pN, uoN are employed to describe unpolarized and polarized data. Bare parameters in 
an effective hadronic Lagrangian are determined in a fit to the data for vrA^ — ► irN, jN — * irN, 
ir~p — ► urn, and 70 — > ujp reactions at center-of-mass energies from threshold to W < 2.0 GeV. The 
T matrix determined in these fits is used to calculate the photon beam asymmetry for w-meson 
production and the loN — ► coN total cross section and ujN scattering lengths. 
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I. INTRODUCTION 



Nucleoli resonances are thought to play a decisive role in reactions of strong, electromag- 
netic and weak probes on nucleons at energies W < 2.0 GeV. The extent to which nucleon 
resonances determine unpolarized and polarized observables in meson production reactions 
and the role of non-resonant contributions in these reactions remains an open question in 
this energy region. Model determinations of the T matrix consistent with the observed 
meson production data in this kinematic regime seeks to resolve the resonance spectrum of 
the nucleon. Such a determination offers the possibility of gaining insight into fundamental 
aspects of quantum chromodynamics, such as the role of chiral symmetry, confinement and 
a detailed understanding of the correlations among the strongly interacting quarks. 

Limited experimental data for to meson production in the resonance region is being rapidly 
augmented. There is existing high precision data from the SAPHIR collaboration [l|] for 
the 7p — > up reaction from which the unpolarized differential cross section (DCS) and 
decay angular distributions have been extracted. Consistent with this data are the more 
recent measurements of the GRAAL collaboration 

as. 

More photoproduction data at 
similar kinematics is anticipated from the CLAS collaboration [4|. The nN — > tuN data 
from bubble and drift chamber experiments is of low precision and there is little overlap in 
different experiment's la, M kinematics. Though there is some discussion in the literature 

u n 

about the validity of the extracted cross sections [3, |8|] we assume the data is correct as 
originally published. 

The importance of including off-shell effects in dynamical coupled channel formulations 
of strong and electromagnetic meson production reactions has been extensively studied 9j, 
[lo| . The present study incorporates off-shell effects in a coupled channel approach and is 
comparable to the model treatments of Krehl et. al. [lOj and Chen et. al. . It should 
be contrasted with coupled channel calculations which take into account coupling of the 
intermediate states only to the continuum and neglect their off-shell contributions such as 
those of the Giessen group [12] and KVI [3]. In the Giessen study, an effective Lagrangian is 
adopted for the channels ttN, inrN, rjN, uN, KA, KH. They assume a resonant contribution 
similar to the one adopted in the present study. 

Motivations for studying the ujN reaction are manifold. Besides the insight into the 
T — i resonance spectrum and implications for meson production reactions, the vector 
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mesons are thought to be important components in very dense matter in the neutron-rich 
stellar environment fl4 . 16]. In nuclear matter, the uNN coupling can play a large role 



in determining the equation of state in some models [17] . 

Any realistic model reaction theory for meson production necessarily incorporates a large 
number of parameters. In the non-resonant terms, we require intermediate state particle 
masses, bare couplings and cutoffs of the hadronic degrees of freedom. For the resonant 
terms we need the bare masses, bare couplings, and cutoffs for each resonance in a given 
partial wave. Given this large number of parameters required to fully determine the T matrix 
the question arises as to the physical and predictive content of such a model. In order to 
address this question in the present study we take the following approach. We determine 
the T matrix by fitting a subset of the available date (described in detail in Section IHip 
and using this, calculate another, unfitted observable (here the photon beam asymmetry 
for u meson production shown in Fig JToT) . The quality of the prediction of the unfitted 
observable is a measure of the utility of model in determining more detailed information 
such as polarization observables from less detailed ones, like unpolarized total or differential 
cross sections. Deficiencies of such a predicted fit indicate the need to improve the model 
assumptions and to augment the included dynamics. 

In the next section we briefly describe the model theory for the six-channel model. The 
results of the fit to the data are presented and discussed in Section IIHI The final section 
gives conclusions and descriptions for improvements to the present study which are under 
development and outlines possible applications of the present approach to other reactions 
including p and production, pion and u meson electroproduction and reactions on nuclear 
targets. 



II. MODEL REACTION THEORY 



The T matrix for ixN — > MB and 7A — * MB - in this work the final MB state is 
restricted to ttN,ujN - is written as a sum of non- resonant, t and resonant, t R contributions 

T(E)=t(E)+t R (E), (1) 

where E = W is the scattering energy of the particles in the center-of-mass frame. Qual- 
itatively, the non-resonant contribution includes rescattering and coupled channel effects 



3 



(a) (b) (c) 

\ _ 

(d) (e) 

FIG. 1: The interactions vm'B',mb and vmb^n include 74 interaction mechanisms. In this figure 
and Figs 121 [3] and H] we show a subset of these. Here the v n pf,irN interaction mechanisms are shown, 
(a) s-channel nucleon exchange; (b) u-channel nucleon exchange; (c) t-channel p exchange; (d) 
s-channel A exchange; (e) u-channel A exchange. 



(a) (b) (c) 

FIG. 2: The VuN,nN interaction mechanisms, (a) s-channel nucleon exchange; (b) u-channel nucleon 
exchange; (c) t-channel p exchange. 

on the Born amplitudes while the resonant contribution includes these effects on the bare 
resonance transition form factors and bare resonance masses. No assumption is made about 
the relative size of the contributions of these terms. Our first objective in this work is to 
determine the T(E) in fits to the observed data. 

Except for the ttN — > ttN reaction where we fit to energy dependent solution of Arndt 
et. al. [18| for the nN — > ixN partial wave amplitudes, we fit the unpolarized and polarized 
cross sections of the 7 and 7r induced reactions. The DCS for it induced reactions is related 
to the T matrix as 

^ = ^PM'*W)p«N{k)\ £ £|TM M , AW ^k',k;£)| 2 , (2) 

M M ,,M N i M N 

where k is the relative momentum of the initial nN state and k' is the relative momentum 
of the final meson-nucleon (M'N f ) state, where M' = n or u. The spin projection of the 
particles in the inital (final) state is {Mm 1 , Mn>). The quantity Pmb(p) = ^EmIpI+EbIp) ' 
where Ei(p) = \/p 2 + mf is related to the density of states. A similar relation holds for the 
unpolarized photoproduction cross section 
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(a) 



(b) 



FIG. 3: The v^n^n interaction mechanisms, (a) s-channel nucleon exchange; (b) tt-channel nucleon 
exchange. 
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TABLE I: The (L, S) terms for partial waves £tj for included channels. 
A. Non-resonant contribution 

The non-resonant contribution to the transition matrix in the partial-wave representation 
for the pion-induced tfT s , M , N , i7rN (E) and the mixed partial- wave/helicity representation for 



(a) (b) (c) 



(d) (e) (/) 

r\/\/\/v 



(5) (ft) (0 

FIG. 4: The v^n^n and v^n^n interaction mechanisms. For photoproduction of ir mesons: (a) 
s-channel nucleon; (b) u-channel nucleon; (c) i-channel ir, p, and u; (d) contact; (e) s-channel A; 
(f) it-channel A. For photoproduction of uj mesons: (g) s-channel; (h) u-channel; (i) i-channel ir 
exchange. 

the photon-induced tjjs'M'N',\ x n t n reactions are 

t J vs'M'N>,e.TTN{k j k; E) = v J I J S i M , N ,^ 1TN {k , k) 

d PP 2tJ J S 'M'N',LSMB( k '^P^ E)G 0>M b(p; E)v J L ^ MB/nN ( P , k), (3) 
^I/S'M'N',\yX N T NlX {^^ Q) = V L?S' M' N' ,\~,\ N T N , Z 0) 

dpp 2 t J Js>M>N>,LSMB( k iP; E)G 0M b(p; E)v j l ^ mb ^ XnTnz (p, q), (4) 

uama - 

where T is the total isospin and T^ jZ is the isospin projection of the nucleon in the initial 
state, J is the total angular momentum of the partial wave, L(L') is the partial wave orbital 
angular momentum of the initial (final) state, S(S') is the total intrinsic spin of the particles 
in the initial(final) state and I = J ± § is the irN initial state orbital angular momentum. 
The included partial waves are shown in Table [fl Channels are defined by the meson species 
M'(M) of the final (intermediate) state and N(N'), the nucleon of the initial (final) state or 
B the intermediate state baryon associated with the meson M of channel MB. The sums, 
J2mb are over channels nN, r]N, it A, <tN, pN, and uN. Here G ,mb(p; E) is the relativistic 
free particle Green's function 

Go ' mb{p] E) = E-E M ( P )-EM-i; MB (p;Ey (5) 

Ei{p) are the free particle energies with masses given in Table [III and Hmb{p', E) is the 
self-energy of the unstable particle in channels MB = nA, aN, pN including the effects 
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FIG. 5: Real part of 7riV — > irN partial wave amplitudes T^ N ^ 7TN (k 7T N , k^N^E) for T = 1/2 versus 
center-of-mass energy, W (in GeV) fit to single energy extraction of Ref. 181 ] ■ The dashed line 
shows the best fit obtained without the second resonance in D15. 

1911. Channels with stable 



induced by the decay of the unstable particle in these channels 
particles MB = nN, t]N, ojN have £a/b = — corresponding to the coupling to the on-shell 
intermediate states. The width of the u meson = 8.5(1) MeV is neglected. 

The vm'B',mb and vmb^n are the effective non- resonant interaction Hamiltonians for 
hadrons tt, r), a, p, lu, N, A, and photons, 7. These interactions are Born amplitudes 
derived from the Lagrangian of Ref.[19| and subjected to the unitary transformation method 
of Ref. {20]. It yields an interaction which is independent of the scattering energy, E = W, 
and depends only on the relative three-momenta of the incoming and outgoing particles. 

In the present model there are 74 interaction mechanisms among the channels 7riV, rjN, 
7rA, aN, pN, uN, and jN. Figures [THU show examples of these explicitly for the terms 
involving 7JV, nN, and uoN channels for pion induced and photoproduction amplitudes. 
In this work we neglect the contribution of the high-mass aj and fj mesons (except the 
/o/er(600)). We make the further simplification in the non- resonant hadronic interaction 
involving the uN channel of including only those terms which couple the uoN to itself and 
to ttN, that is: v^n,mb = v wN ^ n 5 M b,-kN + v w n >u , n 5 M b,u>n- This simplification permits the 
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FIG. 6: Imaginary part of ttN — > ttN partial wave amplitudes T^ N ^N^^N > knN\ E) for T = 1/2 



lfi(| . The dashed 



versus center-of-mass energy, W (in GeV) fit to single energy extraction of Ref. 
line shows the best fit obtained without the second resonance in Di§. 

introduction of a minimal number of additional bare parameters for the u>N channel while 
retaining effects from each of the non-resonant s— , t—, and u— mechanisms. In this way, we 
capture the behavior associated with each mechanism, while maintaining a tractable model. 

Equations and (jlj) represent the bulk of the computational effort required to carry 
out the coupled channel dynamical approach (at the two-body level). Most of the computer 
time required (~ 3/5) is spent evaluating the Born terms. Much of the remainder is spent 
inverting the matrix representing the scattering wave function, T~ x = (1 — vGq)^ 1 = 1 + 



tGo appearing in Eq. 
subtraction methods 



|3| on a momentum grid of 25 Gauss-Legendre points using standard 
2 if ] . Convergence has been checked with grids of up to 45 GL points. 
A parallel FORTRAN90 code has been developed to cope with the long evaluation times for 
a single \ 2 evaluation (~ 100 node-m). It exploits the independence of the partial waves 
and energies in the evaluation of the T matrix. Typically ~ 10 3 — 10 4 y 2 evaluations are 



required for optimization using the minuit package 22 ] . 

The amplitudes tj^ MB tlTN include partial wave contributions up to and including L = 3 
(F-wave). The same is true for the electromagnetic terms except for the t-channel pion 
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J3 

FIG. 7: Real part of 7riV — ► ttN partial wave amplitudes T inN £ nN (k- K N, k^N\E) for T = 3/2 versus 



center-of-mass energy, W (in GeV) fit to single energy extraction of Ref. 18 ] 



exchange in FigJUi). In this case, all partial waves are required for convergence. For L > 3 

the contribution to the electromagnetic non-resonant amplitude tiJs'M>N',\ X n t n z {^i 9! E) 

is calculated at the Born amplitude level only and neglects the effects due to final state 

interactions and coupled channels, ie. the second terms of Eq.(H|. 

The non-resonant interaction depends on the masses of the hadrons and their coupling 

□ 

and cutoff parameters. These values obtained in the five-channel fit of Ref. [191] are shown 
here in Tables HTT1 and HVl for completeness. For the interaction terms, vm'b',mb and vmb^n 
(other than the mass of /o/cr(600) which is a fit parameter), the physical particle masses 
are used. Form factors are included at vertices in the non-resonant interactions, vm'B',mb 
and vmb^n have the form -F(|q|; m) = (A 2 /(A 2 + |q| 2 )) m . Here |q| is either the momentum 
transferred at the vertex or the relative momentum [23]. We use the value m = 2 at all 
vertices. 
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FIG. 8: Imaginary part of ttN — > 7riV partial wave amplitudes T £nN enN (k n N, k n Nj E) for T = 3/2 



versus center-of-mass energy, W (in GeV) fit to single energy extraction of Ref. 
B. Resonant contribution 



18] 



The resonant contribution, t R (E) to the scattering matrix is given as 

B,£irN 

hi 



t?SMBMN(k', k;E) = J2 rfsMB^ (k'; E^iE^^ik; E), ( (5 1 



R,JT 



(k,q;E) = ^ J L T SMB,N : (k;E)V^\E)r J1 



ij y^) 1 - Nj,\-f\ N T N .*Vl- •■ 



(?;£), (7) 



where the sums Sjj run over the resonances in a given partial wave (at most two per channel 
in this work) and Y is the dressed vertex function 

^LSMB,N*(k; E) = Yf^ MBN *(k) 

/ ^ k' 2 t'j^ MB L , SIM , B , (k, k'; E)Go t M'B' {k'; E)Y J y S , M , B , N * (k ) 

L'S'M'B' 



(q;E) = Y N ,^ XNTNz (q) 



+ dk k2T N*,LSMB(k; E)G , M B(k; E)\ 

T O 7V/T JD J 



,JT 
LSMB 



LSMB 



(8) 



x-,\ n t n , z (k, q) (9) 
(10) 
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FIG. 9: Differentia 
data from Refs.[j], 
of each panel. 



cross section vs. center-of-mass angle, 9 for ir p — > ion (in mb/sr) compared to 
. The center-of-mass energy, W (in GeV) is shown in the upper-right corner 



and T>^ (E) is the dressed resonance propagator which depends on the resonance self-energy, 



=:JT 



(11) 

(12) 



V l3 {E) = (E-M^Sy-EyiE) 

X ij (E)= ( k 2 Tj^* LSMB (k] E)Go t MB{k] E)T L ^ SMBN *(k; E). 

LSMB 

The bare vertex functions T should, in principle, be calculated from appropriate ab initio 
hadronic models. This is beyond the scope of the present study. Instead, we parameterize 
the vertex function in the center-of-mass for the partial wave specified by J, L, S for the 
hadronic channels as 

(13) 



^LSMB,N* (fy - (MB 



1 1 r JT k 



In* LSMB 



Here ( M b = -i for MB = ttN, rjN, ttA and ( MB = 1 for M B = aN, pN, uN. At small values 
of the relative MB momentum k, T J L T SMB N , (k) has the form appropriate to the threshold 
production behavior, k L . It is regulated at large k by the form factor, fN*LSMB(k), described 
below. 
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FIG. 10: Unpolarized differential cross section vs. scattering ang 



for 7p — » 7r°p reaction compared to data from Refs. 
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281 ] . The center-of-mass energy, 



W (in GeV) is shown in the upper-right corner of each panel. 



The bare electromagnetic coupling for N* — > 7 A" for all N* except the first P33 resonance 
(number '8' in Tables IVl and IVT1) is given by 

1 



m N 



2tt) 3 / 2 V E N (q R ] 



A JT 

ATV g 




1R JT I \ r 

— 9n*XT n .M> d A,Xv-A 



7— AjV ) 



(14) 



with r$ t _^ t _ XirT 



[-l) J+£+1 / 2 Tf£ X -,x n t Nz where N* is in partial wave £ JT . The 



form for the first P33 resonance is shown in the Appendix. The photon momentum at 



the resonance threshold, is M 



(p) 

N* 



Or 



+ E^{q^) where the resonance mass is taken 



from the Review of Particle Properties'^. The isospin projection of the initial nu- 
cleon is T N ^ Z and the helicities are A 7 and A at. We assume the forms fff*LSMB^) = 
^'n*lsaib I {^n^lsmb + (k — k^*) 2 
kjy* and A^L are fit parameters. 



an d 9 J n*\t n x (?) — 1- The C$ LSMB , A.ff» LSMB , 



III. RESULTS AND DISCUSSION 



The first objective of the present study is the simultaneous description of the pion and 
photon induced single pion and omega meson production data in a coupled channel approach. 
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FIG. 11: Unpolarized differential cross section vs. 9 for ~fp — > 7r + n reaction compared to data from 
Refs. 



25 



26, 



28, 



291 ] . The center-of-mass energy, W (in GeV) is shown in the upper-right corner 



of each panel. 

Recent high precision measurements of uj photoproduction make it possible to strongly 
constrain coupled channel model reaction theories. The DCS and u spin density matrix 
elements (SDME) have been measured at SAPHIR and published 



CLAS collaboration 



1[ and measured by the 



36 1 . The only other observable measured is the single polarization 

nn 

observable, the photon beam asymmetry (PBA) E) at GRA AL [2|, |3| . We have elected 

to also include the older 7r induced reaction data from threshold (~ 1.72 GeV) to 1.764 GeV 

rom the Nimrod synchrotron 6j and the Alvarez detector data from 1.75 GeV to 2.05 GeV 

ij (in 100 MeV bins). 

The world data for pion photoproduction measurements of DCS, shown in Figs fllJlfTTl and 
PBA shown in Figs JT2lfT3l are obtained from the George Washington University Center for 
Nuclear Studies Data Analysis Center [3?}]. These high pr ecision dcitci in 7£> — * tt°p were 
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taken from Refs. 
7P — > 7T + n from Refs. 
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28 1 for DCSandJRe: 
29| ([25 



32 



s. 
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27, 



30 



31 



32 



33[ for PBA and in 



34|) for DCS (PBA), respectively. 



For the purpose of x 2 optimization of the data with respect to the bare parameters of the 
theory a truncated data set of the highest precision data was used which covers as much of 
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FIG. 12: Photon beam asymmetry, Eq(0, E) vs. 9 for ^yp — ► tt°p reaction compared to data from 



Refs. 
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33J. The center-of-mass energy, W (in GeV) is shown in the upper-right 



corner of each panel. 



the angular range as was available on a per energy bin basis. 

In order to accomplish our objective we take as the starting point for this analysis the T 
matrix determined in Ref. 38( which fits the ttN — > ttN partial wave amplitudes 



f JT (k 

J- e«r AT D«r A7 \ ">■ 



(k n N, k n N] E), 



(15) 



extracted from observed data by Ref. 18J in the region 1.1 GeV < E < 2.0 GeV in a five- 
channel approach, excluding uoN . The non-resonant parameters are fixed and shown in 
Tables HT1 IHIl and HVl The resonance parameters (M^l, C$ LSMB , A^ i5MB and k^*) for 
coupling to hadronic channels MB = {iTN,T]N,7TA,aN, pN} shown in Tables IVl and IVT1 



Determination of the six-channel T matrix Tlg MB L , s , M , B ,, T/J mb ^^ AjvTjv proceeds in 
two stages. At the first stage, the ttN — > ttN T = 1/2 partial wave amplitudes of Figsj5] 
and El the n~p — > um DCS of Figj9] and the 7p — > up DCS of FigJTU are fit simultane- 



t 

uNN 



ously. This is accomplished by adjusting the non- resonant couplings glj NN , ^nn an d A 
appearing in the s— and u— channel u> emission and absorption of Figsj2]^a),(b), Figsj3]and 
FigsjU(g),(h) and by adjusting the resonance parameters N* — > uN, Gj}, LSujN and A.fJ* LStjN - 
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3eam asymmetry, E) vs. 9 for jp — ► 7r + n reaction compared to data from 



32, 



33 



341 ] . The center-of-mass energy, W (in GeV) is shown in the upper-right 



corner of each panel. 



The introduction of the ujN channel to the calculation requires the addition of a second -D 15 
resonance, shown in bold type in Tables [V] and IVIt in order to fit the data. These points 
will be discussed in more detail below. 

At the second stage of the fit, all non- resonant and hadronic channel resonant parameters 
are fixed and the single meson photoproduction data is fitted. Pion photoproduction data 
used for the fit includes the DCS in Figs{10] and [11] and the PBA in Figs{T2] and [13] in 
the region 1.1 GeV < E < 2.0 GeV. Omega meson photoproduction data used for the fit 
includes only the SAPHIR measurement [l| of the DCS from threshold, 1.72 GeV to 2.0 GeV 
shown in FigJTBI This is accomplished by varying the photon helicity couplings, A^, XTn for 
A = |, § and T NjZ — +i. The resulting fits compared to the existing world data are shown 
as solid curves in Figs fTUHTBl 

The overall quality of the fits to the complete set of data are in fair agreement for energies 
E < 1.65 GeV. The T = 1/2 ttN — > ttN partial wave amplitudes of Figj5] agree at the 1 — a 
level for all partial waves except the two highest. The T = 3/2 partial fits are of similar 
quality except for the S\\ wave and the P\\ wave at energies E > 1.9 GeV. The fits to the 
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FIG. 14: Unpolarized differential cross section for ~fp — > cop versus 9 (in //b/sr) compared to data 
from Ref.[l|. The center-or-mass energy, W (in GeV) is shown in upper-right of each panel. 



photoproduction data are good at low energies but degrade significantly at E > 1.65 GeV 
especially in the 7p — > n°p reaction. Coupling to the inrN channel is expected to be large 
here. 

Large values of g U NN have been deduced from studies of the nucleon electromagnetic form 



factors 



391 ] and various NN studies 



These studies yield a range of 10 < g U NN % 



42] though with a strongly 



20. Values of this order were assumed for the studies in Ref. 
suppressing form factor due to a small cutoff, Aotl = 0.5 GeV. At early stages of the fit 
when we attempted to use the values 5^ aw, k Uj and A^nn determined in fits to the nN — > nN 
data the resulting cross sections were too large by one or two orders of magnitude for both 
7i~p — > um and 7p — > up reactions. In order to reproduce the data within the present model 
for the limited parameter search which we have performed it was necessary to introduce the 
non- resonant coupling parameters g^ NN , k^, and A^^y. These parameters appear at vertices 
in the graphs of Figsj2]^a),(b), FigsJH] and Figs]4|g),(h). The interactions correspond, in a 
four- dimensional formulation to vertices with timelike momentum transfer. In the fits to 
the ttN — ► ttN data uj mesons appear only in graphs corresponding to spacelike momentum 



transfer. An similar situation obtains in pp — > ppn° reactions 43j where different couplings 
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FIG. 15: Predicted photon beam asymmetry, Y, w {0, E) for jp — > top (solid curve) for E = W shown 
in lower- left corner of each panel compared with data from Ref . [2] . At the lowest energy, the effect 
of removal of various resonances is shown. Removing all ('nr only' - thin dashed) and Su (dotted), 
.D13 (dashed) , F15 (dot-dashed). 



are used for exchanged and emitted s-wave pions. The small value obtained for the g^^N 
coupling is near the result found by the Giessen group's study [12]. The treatment here 
is certainly phenomenological but no more so than introducing other non-resonant reaction 
mechanisms involving heavy mesons (eg. including v u jf,pN or f /a exchange in Vun,*n) or 
resonances. These alternatives should nevertheless be explored as a guide, at least, to the 
model dependencies of the present approach. 

We have found that the introduction of the usN channel significantly modifies the behavior 
of the tcN —+ ttN _D 15 partial wave amplitude. This can be seen in Figs|S] and M where we 
show as a dashed-line curve in the -D15 panel the optimal curves found in the first stage fit 

to nN — > tyN, ix~p — > un, and jp — > up data described above. The A 5 ^ photocoupling 

2 

is large and could be an important effect in, for example, the electroproduction reaction. 



Comparison of the iV* — > ujN physica 



masses and branching fractions determined in this 



work with other calculations (as in Ref. 35j) require the analytic continuation of the T matrix 



amplitudes to the physical pole position; this work is in preparation. 
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m,N 938.5 
138.5 
547.5 



m A 1300.0 
898.6 
811.7 
782.6 



m 



in 



m 



TABLE II: Propagator masses 
(MeV) appearing in Eq.([5]). 



fnNA 
fr)NN 
9pNN 
Kp 

9u>NN 

K L 

9oNN 

9pTT1T 

fwAA 
fpNA 

9c7V7V 
9ulTTp 

9pAA 
KpAA 

9piT*f 
9iOTT'y 

m n 



0.08 
2.206 
3.889 
8.721 
2.654 
8.100 
1.020 
1.298 
1.002 
6.815 
4.000 
1.000 
7.516 
2.353 
6.956 
3.302 
2.000 
0.1027e 
0.3247e 
500.1 MeV 





810 


A-ttNA 


829 




1087 


Ap7T7T 


1094 




1523 




589 




624 




781 


A-pNA 


1200 


Attaa 


600 


A,j7r7r 


1200 


Aj7rp 


600 


ApAA 


600 



TABLE IV: Lagrangian non- 
resonant cutoffs (MeV). 



TABLE III: Lagrangian bare 
coupling and a mass. 



The prediction for the PBA [44J in 7p — > up 



E) 



(16) 



is shown in Figjl~5l Here crii(cr_|_) is the differential cross section for linearly polarized photons 
in (perpendicular to) the emission plane of the u meson. At the lowest energy E = 1.743 
GeV a study is made of the sensitivity to the resonance contribution for three cases. The 
thin-dashed line is the result when all the resonance contributions have been removed. Other 
curves in the figure show the result when one of three dominate waves is removed. 
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# L TJ 


M(°) 




nN 


7/iY 


7rA 




aN 




pN 










Ai 


A 3 

■> 


1 


Sxx 


1800 


99.9 


7.049 


9.100 


-1.853 




-2.795 


2.028 


0.027 




-3.761 


0.405 




83.8 




2 


Sxx 


1880 


100.0 


9.824 


0.600 


0.045 




1.139 


-9.518 


-3.014 




-0.516 


0.366 




-40.3 




3 


Szx 


1850 


20.7 


5.275 




-6.175 






-4.299 


5.638 










129.4 




1 


Pxx 


1763 


76.1 


3.912 


2.621 


-9.905 




-7.162 


-5.157 


3.456 




-3.362 


5.231 




-21.8 




5 


Pxx 


2037 


22.1 


9.998 


3.661 


-6.952 




8.629 


-2.955 


-0.945 




-2.095 


1.043 




-27.5 




6 


PX3 


1711 


76.4 


3.270 


-0.999 


-9.988 - 


5.038 


1.015 


-0.003 


2.000 


-0.081 


5.737 


-0.548 


-0.204 


-12.4 


-63.8 


7 


P3X 


1900 


100.0 


6.803 




2.118 






9.915 


0.153 










54.1 




8 


P 33 


1603 


83.9 


1.312 




1.078 


1.524 




2.012 


-1.249 


0.379 








-78.6 


-131.2 


9 


P33 


1391 


-93.3 


1.319 




2.037 


9.538 




-0.317 


1.036 


0.766 








-6.7 


5.3 


10 


D ri 


1899 


-35.3 


0.445 


-0.017 


-1.950 


0.978 


-0.482 


1.133 


-0.314 


0.179 


-0.081 


3.740 


0.230 


88.8 


-71.4 


11 


D 13 


1988 


-41.7 


0.465 


0.357 


9.919 


3.876 


-5.499 


0.289 


9.628 


-0.141 


7.883 


9.900 


3.386 


-54.5 


46.8 


12 


Dxs 


1898 


0.0 


0.312 


-0.096 


4.792 


0.020 


-0.455 


-0.179 


1.249 


-0.101 


0.625 


1.086 


-0.156 


33.0 


40.3 


13 


D15 


2334 


9.7 


0.167 


-0.106 


0.190 - 


0.098 


-0.075 


-0.530 


0.228 


0.099 


-0.150 


-1.990 


0.199 


12.6 


87.4 


11 


D33 


1976 


36.7 


0.945 




3.999 


3.997 




0.162 


3.949 


-0.856 








95.9 


-6.1 


15 


Fis 


2187 


92.1 


0.062 


0.000 


1.040 


0.005 


1.527 


-1.035 


1.607 


-0.026 


-0.046 


2.212 


0.078 


-99.8 


-68.1 


16 


F35 


2162 


-84.2 


0.174 




-2.961 - 


1.093 




-0.076 


8.034 


-0.061 








-61.0 


-103.4 


17 


F37 


2137 


-100.0 


0.254 




-0.316 - 


0.023 




0.100 


0.100 


0.100 








45.9 


47.7 



TABLE V: Bare masses Af<°) (MeV) appearing in the resonance propagator of Eq.([5]), and 
the ranges k^* (MeV), strong couplings Gj^ UB N * (MeV~ 1//2 ) and photo-couplings Aj^ (1CT 3 
GeV" 1 / 2 ) in Eqs.CCl and CE}. 



# 




nN 


TjN 




7rA 




aN 






pN 








luN 




1 


Sxx 


1676 


4 


599.0 


554.0 






801.0 


1999 


9 


1893.7 




500 


2 


817.9 




2 


Sxx 


533 


5 


500.0 


1999.1 






1849.5 


796 


8 


500.0 




503 


1 


622.0 




3 


S31 


2000 







500.0 








500 





500.0 












■1 


Pxx 


1203 


6 


1654.8 


729.0 






1793.1 


622 





1698.9 




675 


8 


516.9 




5 


Pxx 


646 


9 


897.8 


501.3 






1161.2 


500 


1 


922.3 




533 


7 


950.1 




6 


PX3 


1374 





500.2 


500.0 


500 


8 


640.5 


500 





500.1 


1645.2 


500 


1 


547.3 


513.3 


7 


Psx 


828 


8 




2000.0 








1998 


8 


2000.0 












8 


P33 


746 


2 




846.4 


781 







585 





500.2 


1369.7 










9 


P33 


880 


7 




507.3 


501 


7 




606 


8 


1043.4 


528.4 










10 


Dxs 


1658 





1918.2 


976.4 


1034 


5 


1315.8 


599 


8 


1615.1 


1499.5 


565 


5 


802.9 


978.1 


11 


D 13 


1094 





678.4 


1960.0 


660 





1317.0 


550 


1 


597.6 


1408.7 


500 


5 


506.2 


545.2 


12 


D 15 


1584 


7 


1554.0 


500.8 


820 


2 


507.1 


735 


4 


749.4 


937.5 


1036 





999.2 


996.0 


13 


D15 


1223 


8 


1990.2 


1910.4 


996 


1 


921.6 


1022 





1941.9 


997.0 


930 


2 


998.2 


999.2 


11 


D33 


806 







1359.4 


608 


1 




1515 





1999.0 


956.6 










15 


Fxs 


1641 


6 


655.9 


1899.5 


522 


7 


500.9 


500 


8 


500.0 


1060.9 


541 


8 


502.0 


651.6 


16 


F35 


1035 


3 




1228.0 


586 


8 




1514 


8 


593.8 


1506.0 










17 


F37 


1049 







1180.2 


1031 


8 




600 





600.0 


600.0 











TABLE VI: Resonance strong form factor cutoff parameters in MeV. 
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The total cross section of the reaction uN — > uN is of interest for realistic calculations 
of nuclear matter properties. The predicted total cross section for this reaction is shown in 
FigJTBl The scattering lengths obtained from the T matrix are 



aj 



lim 



vrm^mjv j 



ai 

2 



II > 

2 



(17) 

(18) 
(19) 



[-0.0454 - i0.0695] fm, 
[0.180 -i0.0597] fm, 

related to the total cross section at threshold by a^N^E — > m^ + m^v) = 47r(|ai | 2 + 2|a|| 2 )/3. 

The total cross section for ^p — > up is shown in FigJTT] along with contributions from 
partial waves irj with significant contributions. The error bars on the total cross section 
are statistical jlj . Systematic errors are shown in Ref.pJ to be about 10 — 15%. They arise 
from, among other sources, the extrapolation of the DCS in the forward and backward di- 
rections for center-of-mass a;-meson scattering angles 9 < 15° and 9 > 150°. The systematic 
errors for the DCS from Ref . [1] are largest in the backward direction where the discrepancy 
from our calculated cross section, as seen in Fig{18] is most pronounced. Nevertheless, the 
present calculated — > up DCS appears to miss some small angle structure near x ~ 0.5 
above E > 1.86 GeV and x ~ —0.5 and 0.5 for E > 1.935 GeV, possibly a result of de- 
structive interference effects of some combination of additional higher mass resonances and 




FIG. 16: Calculated total elastic tuN cross section (mb) as a function of center-of-mass energy 
E = W. 
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FIG. 17: Total cross section for — ► up (/ib) compared with that extracted from data from 



Ref. 



451 ] (circles) and Ref. [1(1 (squares) as a function of center-of-mass energy E = W . The error 
bars on the data from Ref.pj are statistical errors only. Systematic errors are about 10 — 15% Qj. 
The three partial waves with the largest contribution are shown. See text for discussion. 



non-resonant effects not included in this study. 49] 



IV. CONCLUSION 



A dynamical coupled channel model for six-channels has been employed in simultaneous 
fits of the pion and photon induced single pion and omega production reactions. The ttN 
partial wave amplitudes, unpolarized differential cross sections (DCS) and photon beam 
asymmetry (PBA) have been fit with % 2 ~ 1 for center-of-mass energies from threshold 
to E < 1.65 GeV. At higher energies, the model is unable to accurately reproduce the 
data. There are several ways one might attempt to remedy this deficiency. If we work 
within the present model formulation and keep the same non-resonant mechanisms it is 
possible that a more thorough search of the parameters may yield a better fit at higher 
energies. The introduction of more resonances may also yield a better fit. However, for 
energies W > 1.8 GeV the existing data makes distinguishing non-resonant mechanisms 
from resonant mechanisms difficult. Particularly useful in this endeavor would be more 
high precision single and some double polarization observables for both pion and cu-meson 
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01 I I o.i I ' l I 

-1 -0.5 0.5 1 -1 -0.5 0.5 1 



COS 

FIG. 18: Semi-logarithmic version of Figll44 unpolarized DCS for jp —* ujp (^b/sr) vs. cos#. 
Center-of-mass energies are shown in the upper-right corner of each panel. The sharp angular 
features near are accenuated x ~ 0.5 for E > 1.86 GeV and x — —0.5 for E > 1.935. 

production. 

Additional mechanisms in the non-resonant terms are sure to contribute, perhaps signif- 
icantly, to the calculated scattering observables at these higher energies. At the two-body 
level we have neglected couplings such as v^n^n, v^n^a, v U N,aN, an d v u ^ )P n- There may 
also be significant effects from additional mechanisms in the i^jv^jv interaction. We have 
neglected the effects of t-channel a exchange, T] exchange (generally thought to be small) and 
Pomeron exchange in 7p — > up, known to have large contributions at forward angles at high 
energies. At energies above the two-pion production threshold, the ttttN channel contribu- 
tion can give a significant contribution and must be calculated. This can be accomplished 
in the present model formulation and is currently under study. 

A fair prediction for the u meson PBA, E) has been obtained near threshold, shown 
in the upper-left panel in FigJT5l At higher energies, W > 1.8 GeV the calculated beam 
asymmetry does not agree well with the shape of the data. This observable is sensitive to 
both resonant and non-resonant contributions could improve with any of the refinements 
discussed above. 
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The present model will be used to analyze the electroproduction data in the region Q 2 < 5 
GeV 2 and the extension to the photoproduction of p and </> vector mesons. Photoproduction 
of to from nuclei gives information about the spectral function of the uj meson in the nuclear 

n 

medium[46j. Inclusion of off-shell effects in such an analysis is required and the present 
model affords a starting point for their inclusion. 



APPENDIX: P 33 (1232) TRANSITION FORM FACTOR 



The transition form factor for 7^ — > A (1232) is taken as 47] 



Tf(q) = -K. M G M (q 2 ) + (/C M + lC E )G E {q 2 ) 

2 >P 

Tf p (q) = * [-£mGm(? 2 ) + (/Cm - JC E )G E (q 2 )} 



/Cm 
JC e 



'E N (q)+m N 1 3(m A + m N ) E\q\<fK 



(2vr)3/2y 2E N (q) ^2\q^\ 2m N Q 2 + (m A + m N ) 
4E\q\ 2 1 



r/C 



M ■ 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



E N (q) + m N Q 2 + (m A - m N ) 2 
On resonance at the photon point, Q 2 = the Gm(0) and G E (0) are related to the value 
Af n in Table |V] as 

y/3 i /mA|q| 



A 33 



hp 



Af 

3.P 



2 2rriN \] 

1 e /?7T,A|q| 



[G M (0) + G E (0)} 



2 2mjy V n^N 
The values in the table correspond to 



[G M (0)-3G E 



Gm(0) = 1.62 
G B (0) = 0.015. 



(A.5) 
(A.6) 



(A.7) 
(A.8) 
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